Resolved X-ray emission line profiles 
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By quantitatively fitting simple emission line profile models that include both atomic opacity 
and porosity to the Chandra X-ray spectrum of C Pup, we are able to explore the trade-offs 
between reduced mass-loss rates and wind porosity. We find that reducing the mass-loss rate 
of £ Pup by roughly a factor of four, to 1.5 x 10 -6 M Q yr _1 , enables simple non-porous wind 
models to provide good fits to the data. If, on the other hand, we take the literature mass-loss 
rate of 6 x 10 -6 M Q yr _1 , then to produce X-ray line profiles that fit the data, extreme porosity 
lengths - of w 3 R, - are required. Moreover, these porous models do not provide better fits 
to the data than the non-porous, low optical depth models. Additionally, such huge porosity 
lengths do not seem realistic in light of 2-D numerical simulations of the wind instability. 



1 Introduction 

The X-ray emission lines from stars are wind- 
broadened but are surprisingly symmetric. Asym- 
metry should arise due to the preferential absorp- 
tion of red-shifted photons, which are produced on 
the far side of the wind and seen through a larger col- 
umn density of cold, absorbing material. The very 
modest asymmetry in the observations implies that 
mass-loss rates are lower than has been presumed, 
which is, of course, in line with other recent observa- 
tions. Alternatively, it has been suggested that the 
mass-loss rates are actually high, but the effective 
opacity of the wind is reduced due to porosity, or 
macroclumping (Oskinova et al. 2006). 

In this very brief paper, in which we focus on one 
strong representative line in the Chandra grating 
spectrum of £ Pup, we explore whether the mea- 
sured profile shapes can discriminate between mass- 
loss reduction and porosity. Even if they cannot, by 
fitting models to data we can quantitatively explore 
the trade-off between the key parameters: fiducial 
optical depth, r* = kM /47rz; 00 i? H », and the terminal 
porosity length, hoo- Here we do this by fitting the 
Owocki & Cohen (2001) profile model to one strong 
line in the Chandra HETGS/MEG spectrum of ( 
Pup. We then fit a modified profile model, where the 
opacity is adjusted for porosity according to Owocki 
& Cohen (2006). 



2 Fitting Models to the Data 

There are 560 counts in the Fe XVII line shown in 
Fig. 1, accumulated during an exposure time of 68 
ks. The line is well resolved and clearly asymmetric, 



with the characteristic blue shifted and skewed line 
shape that is expected from a spherically symmetric, 
beta-velocity wind with embedded X-ray emission 
and continuum absorption from the cool, dominant 
component of the wind. The first model we fit to the 
data is the simplest - a non-porous Owocki & Co- 
hen (2001) model where the optical depth depends 
on the atomic opacity. We note that this model does 
not preclude microclumping of the sort that affects 
density-squared emission diagnostics, but does not 
affect column-density based diagnostics like X-ray 
emission line profiles. However, this standard pro- 
file model does not include large scale porosity, or 
macroclumping, with h = £/ f > r, where h is the 
porosity length, defined as the clump size scale to 
the volume filling factor (Owocki & Cohen 2006). 
We do fit porous models later in this paper. 

We fit the model, along with a fixed power-law 
continuum, within XSPEC, to the line over the 
wavelength range shown in Fig. 1. The fitting rou- 
tine in XSPEC adjusts the free parameters of the 
model: i?o, the radius below which there is no X-ray 
emission, r* , the fiducial wind optical depth, and the 
normalization, until the global minimum of the fit 
statistic - here the C statistic (Cash 1979), which 
is the maximum likelihood statistic for Poisson- 
distributed data - is found. The best-fit model, with 
i?o = 1.5 ± .2 R* and t* = 2.0 ± .5, is plotted over 
the data in Fig. 1. 

The uncertainties on the parameters can be esti- 
mated by evaluating AC, the difference between the 
C statistic value for a given choice of parameters and 
that of the best-fit model. The distribution of AC 
is the same as that of A% 2 and a specific value cor- 
responds to a formal confidence limit. By drawing 
contours of constant AC in the parameter space of 
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interest, we can see the extent of a given parameter's 
possible range and the correlation between parame- 
ters, as we show in Fig. 2. 
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Figure 1: An Fe XVII emission line, formed in a 
thermal plasma at a few 10 6 K, with Pois- 
son error bars. The approximate spectral 
resolution is indicated by the orange bar. 
The best-fit non-porous model is shown. 



M = 1.5±.4xl0 6 M yr 1 , given an X-ray opacity 
at 15 A of 70 cm 2 g" 1 (Waldron et al. 1998). If, on 
the other hand, we take the mass-loss rate of ( Pup 
to be M = 6 x 1(T 6 M yr" 1 and use R* = 19 R Q 
and Voo — 2250 km s _1 (Puis et al. 1996), then we 
expect the wind to have a fiducial optical depth of 
t* = 8. The best- fit model is favored over ar, = 8 
model with > 99.99% confidence (AC = 79). In Fig. 
3 we show the optimal model with a fixed t* = 8 (i?o 
and normalization are free parameters) along with 
the global best-fit model. 
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Figure 2: These confidence contours in the Rq - t* 
parameter space of the fit of the non- 
porous model to the Fe XVII line, shown 
in Fig. 1, enclose 68%, 90%, and 95% of 
the model parameter space where the true 
model parameters lie, given the data and 
the assumed model. 



The value of r* 
ting, r* = 2.0 ± 



we derive from the model fit- 
5, implies a mass-loss rate of 



Figure 3: The global best-fit non-porous wind model 
is shown in red and the best-fit model with 
the wind opacity expected based on the 
standard mass-loss rate is shown in blue. 




Figure 4: Suite of porous line profile models, shown 
at infinite resolution. Within each panel, 
three optical depth values are represented, 
t* = 1,2,8: black, red, blue. 
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Wc next fit models that allow for porosity, but are 
otherwise identical to the Owocki & Cohen (2001) 
models fit in the previous subsection. In these 
porous models, the opacity is not atomic, but rather 
geometric, and due to optically thick clumps of size 
i and filling factor /. The optical depth is thus mod- 
ified according to K e ff = K (l — e ~ Tc )/ r e, where k is 
the microscopic opacity and r c is the clump optical 
depth, t c = n£(p)/f, where (p) is the local mean 
wind density (Owocki & Cohen 2006). We parame- 
terize the porosity length as h(r) = h^l — R*/r)P 
and fix (3 = 1, leaving the terminal porosity length, 
/ioo as the single additional parameter of the model. 
As can be seen in Fig. 4, porosity has only a weak 
effect, even at hoo = 1. Values above unity are re- 
quired before the profiles are strongly affected. 

Now, by fitting this porous line profile model to 
the data, we can examine the joint constraints on the 
key parameters: r» and hoo, while allowing Rq and 
the normalization to also be free parameters of the 
model fit. When we do this, we find a best-fit termi- 
nal porosity length of hoo = 0.0. In other words, the 
non-porous model is preferred over the porous one. 
In Fig. 5 we show the joint confidence limits on hoo 
and r». 




Figure 5: The 68%, 90%, and 95% confidence con- 
tours in t* - hoo parameter space, with the 
best-fit global model shown as the red star. 
The green star is the location of the opti- 
mized model with r* = 8 fixed at the value 
implied by the literature mass-loss rate. 
The green line denotes the 68% confidence 
range on the terminal porosity length, hoo 
when t* = 8. The vertical red line em- 
phasizes that even with = 1, the wind 
optical depth is only a little higher than 
the non-porous value of t* = 2.0. 



From this figure, we can see that the model with 
low optical depth and no porosity is preferred with 
a high degree of formal significance over the model 
with high optical depth and a large porosity length. 
We can also see, from the shape of the confidence 
contours, that the porosity length does not apprecia- 
bly change the optical depth derived from the data 
until hoo > R*- This is in accord with the expec- 
tations from the sensitivity of the model profiles to 
the value of this parameter (as shown in Fig. 4). 

We stress that the porous model with the higher 
mass-loss rate provides an adequate fit to the data, 
just not as good a fit as the non-porous model. It 
is possible that a porous model with a different de- 
scription of porosity could fit the data as well as 
the non-porous model fit we show. However, we can 
definitively say two things from the model fitting we 
have reported on here: (1) Porosity cannot explain 
the only modestly asymmetric line profiles in £ Pup 
any better than a non-porous model with a reduced 
mass-loss rate can; if anything, it provides a worse 
fit. And, (2) in order to fit the data with a high 
mass-loss rate, high porosity model, terminal poros- 
ity lengths of several R* are required. 

Regarding this last point, two-dimensional sim- 
ulations of the line-driven instability show numer- 
ous very small structures, as the compressed shells 
created by the instability break up laterally. Con- 
straints on moving emission bumps seen in WR 
winds using these 2-D simulations shows that the 
clumps subtend about 3 degrees, or 1/20 steradian, 
implying linear clump scales of 1/20 the local stel- 
lar radius, and thus quite small porosity lengths 
(Dessart & Owocki 2003, 2005). 
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